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ABSTRACT 

Aiming at exploring the nature of dark energy, we use thirty-six observational Hubble parameter data (OHD) 
in the redshift range 0 < z < 2.36 to make a cosmological model-independent test of the two-point Omh^(z 2 ', Zi) 
diagnostic. In ACDM, we have Omh^ = where is the matter density parameter at present. We bin 

all the OHD into four data points to mitigate the observational contaminations. By comparing with the value 
of which is constrained tightly by the Planck observations, our results show that in all six testing pairs of 
Omh^ there are two testing pairs are consistent with ACDM at Icr confidence level (CL), whereas for another 
two of them ACDM can only be accommodated at 2cr CL. Particularly, for remaining two pairs, ACDM is 
not compatible even at 2cr CL. Therefore it is reasonable that although deviations from ACDM exist for some 
pairs, cautiously, we cannot rule out the validity of ACDM. We further apply two methods to derive the value of 
Hubble constant //q utilizing the two-point (9m/!^(z2;zi) diagnostic. Weobtain//o = 71.23+1.54 km s“' Mpc“* 
from inverse variance weighted Omir value (method (I)) and Hq — 69.37 + 1.59 km s * Mpc * that the Omh^ 
value originates from Planck measurement (method (II)), both at Icr CL. Einally, we explore how the error 
in OHD propagate into w(z) at certain redshift during the reconstruction of w(z). We argue that the current 
precision on OHD is not sufficient small to ensure the reconstruction of w(z) in an acceptable error range, 
especially at the low redshift. 

Subject headings: cosmological parameters — Hubble constant — Hubble parameter: error — dark energy: 
error — methods: statistical 


1. INTRODUCTION 

In the past few decades, a number of approaches have been 
well established to quantitatively study the expansion history 
and structure growth of the universe (see Erieman et al. 2008; 
Weinberg et al. 2013, for recent reviews). Originally aimed 
at measuring the cosmic deceleration rate (Riess et al. 1998; 
Perlmutter et al. 1999), the observations of Type la super¬ 
novae (SNIa) are instead providing ample evidence for an ac¬ 
celerating expansion to an increasing precision (Suzuki et al. 
2012). Tliis acceleration is also strongly supported by other 
complementary probes, including the measurements of the 
baryon acoustic oscillation (BAO) features (Eisenstein et al. 
2005), the weak gravitational lensing (Kilbinger et al. 2013), 
the abundance of galaxy clusters (Benson et al. 2013), the 
cosmic microwave background (CMB) anisotropies (Planck 
Collaboration et al. 2015), the linear growth of large-scale 
structure (Dodelson et al. 2002), and the Hubble constant Hq 
(Ereedman et al. 2012). As the primary motive of modern 
cosmology shifts from whether the universe is accelerating to 
why, a robust physical model that explains the cosmic accel¬ 
eration is still under debate. 

A prevailing interpretation is the existence of an exotic en¬ 
ergy constituent - often coined dark energy (DE) - with neg¬ 
ative equation of state (EOS) parameter w = Pde/Pde- By 
far, the most popular model for DE remains to be the sim¬ 
ple cosmological constant cold dark matter model (ACDM), 
with w = -1 at all cosmic periods (Li et al. 2011). However 
as popular as the ACDM model is, it still suffers from the 
fine tuning and coincidence problems (Weinberg 1989; Zlatev 
et al. 1999). In addition, it has been noticeably argued about 
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the possibility of DE evolving its EOS, i.e. the dynamical DE 
models (w = w(z)), amongst which there exit Quintessence 
(w > -1, Ratra cfe Peebles 1988; Wetterich 1988), Phantom 
(w < -1, Caldwell 2002), K-essence (w > -1 orw < -1, 
Armendariz-Picon et al. 2000, 2001), and especially Quintom 
(w crossing -1, Peng et al. 2005, 2006) models. Nonetheless, 
all these models still await more physically motivated, pro¬ 
found understanding. In the meantime, it is crucial to estab¬ 
lish tests which are based upon direct observations and ca¬ 
pable of unveiling, if any, dynamical features of DE. One of 
these diagnostics is Om{z), which is defined as a function of 
redshift z (Sahni et al. 2008; Zunckel & Clarkson 2008), i.e., 


Om(z) - 


h\z) - 1 
(l+z)3-l’ 


( 1 ) 


with h - and H(z) denoting the Hubble expansion rate. 
Om(z) has a property of being - the matter density param¬ 
eter at present - for oj — -1. Moreover, Shafieloo et al. (2012) 
revised this diagnostic to accommodate two-point situations, 
i.e.. 


Om(z2;zi) 


P(Z2) - h\zi) 

(1 +Z2)^ - (1 +Zl)^' 


( 2 ) 


In this case, the measurement of Om(z2', Zi) = Om is a remark¬ 
able piece of evidence for an underlying ACDM model. In 
other words, the measurement of Om{z2',Zi) + Dm implies a 
deviation from ACDM and the fact that an evolving EOS with 
redshift should be considered. Consequently, according to the 
observations of the cosmic expansion history, we can success¬ 
fully distinguish between DE models and perform correlative 
studies in cosmology using this diagnostic. 

In this paper, we first summarize the observational methods 
for Hubble parameter measurements and introduce the cur¬ 
rently available data sets in Section 2. In Section 3, using the 
observational H(z) data (OHD) hereinafter mentioned, we test 
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the ACDM model, according to the diagnostic defined above. 
We also use this diagnostic to derive Hq values in this section. 
In Section 4, we turn to trace out the effect of OHD errors on 
the reconstruction of DE EOS. Einally, we summarize our re¬ 
sults and discuss limitations as well as prospects in Section 5. 

2. THE OBSERVATIONAL H(Z) DATA SETS 

OHD can be used to constrain cosmological parameters be¬ 
cause they are obtained from model-independent direct ob¬ 
servations. Until now, two methods have been developed to 
measure OHD: galaxy differential age and radial BAO size 
methods (Zhang et al. 2010). Jimenez & Loeb (2002) first 
proposed that relative galaxy ages can be used to obtain H{z) 
values and they reported one H{z) measurement at z ~ 0.1 
in their later work (Jimenez et al. 2003). Simon et al. (2005) 
added additional eight H{z) points in the redshift range be¬ 
tween 0.17 and 1.75 from differential ages of passively evolv¬ 
ing galaxies, and further constrained the redshift dependence 
of the DE potential by reconstructing it as a function of red¬ 
shift. Stern et al. (2010) provided two new determinations 
from red-envelope galaxies and then constrained cosmologi¬ 
cal parameters including curvature through the joint analysis 
of CMB data. Eurthermore, Moresco et al. (2012) obtained 
eight new measurements of H{z) from the differential spectro¬ 
scopic evolution of early-type, massive, red elliptical galaxies 
which can be used as standard cosmic chronometers. Using 
luminous red galaxies from Sloan Digital Sky Survey (SDSS) 
Data Release 7 (DR7), Chuang & Wang (2012) measured a 
new H{z) point at z = 0.35. Later, by applying the galaxy 
differential age method to SDSS DR7, Zhang et al. (2014) 
expanded the H{z) data sample by four new points. Tak¬ 
ing advantage of near-infrared spectroscopy of high redshift 
galaxies, Moresco (2015) obtained two latest measurements 
of//(z). 

In addition, H{z) can also be extracted from the detection 
of radial BAO features. Gaztanaga et al. (2009) first obtained 
two H{z) data points using the BAO peak position as a stan¬ 
dard ruler in the radial direction. Blake et al. (2012) further 
combined the measurements of BAO peaks and the Alcock- 
Paczynski distortion to find three other H{z) results. Samushia 
et al. (2013) provided a H{z) point at z = 0.57 from the BOSS 
DR9 CMASS sample. Xu et al. (2013) used the BAO signals 
from the SDSS DR7 luminous red galaxy sample to derive 
another observational H{z) measurement. The H{z) values 
determined based upon BAO features in the Lyman-aforest 
of SDSS-III quasars were presented by Busca et al. (2013), 
Eont-Ribera et al. (2014), and Delubac et al. (2015), which 
are the farthest observed H{z) results so far. The above data 
points are all presented in Table 1 and marked in Eigure 1. 

We use a ACDM model with no curvature term to compare 
theoretical values of Hubble parameter with the OHD results, 
with the Hubble parameter given by 

H(z) = Hq Vilm(l+z)3 + l-0„, (3) 

where cosmological parameters take values from the Planck 
temperature power spectrum measurements Planck Collab¬ 
oration et al. (2014). The best fit value of Hq is 67.11 
km s ' Mpc *, and is 0.3175. The theoretical computa¬ 
tion of H{z) based upon this ACDM is also shown in Eigure 1. 

Being independent observational data, H{z) determinations 
have been frequently used in cosmological research. One of 
the leading purposes is using them to constrain DE. Jimenez 
& Loeb (2002) first proposed that H{z) measurements can be 


TABLE 1 

The currently available OHD measurements 


z 

mz) 

Method 

Ref. 

0.0708 

69.0 ± 19.68 

I 

Zhang et al. (2014) 

0.09 

69.0 ± 12.0 

I 

Jimenez et al. (2003) 

0.12 

68.6 ± 26.2 

I 

Zhang et al. (2014) 

0.17 

83.0 ±8.0 

I 

Simon et al. (2005) 

0.179 

75.0 ±4.0 

I 

Moresco et al. (2012) 

0.199 

75.0 ±5.0 

I 

Moresco et al. (2012) 

0.2 

72.9 ± 29.6 

I 

Zhang et al. (2014) 

0.240 

79.69 ± 2.65 

II 

Gaztanaga et al. (2009) 

0.27 

77.0 ± 14.0 

I 

Simon et al. (2005) 

0.28 

88.8 ± 36.6 

I 

Zhang et al. (2014) 

0.35 

82.Uj| 

I 

Chuang et al. (2012) 

0.35 

84.4 ± 7.0 

II 

Xu etal. (2013) 

0.352 

83.0 ± 14.0 

I 

Moresco et al. (2012) 

0.4 

95 ± 17.0 

I 

Simon et al. (2005) 

0.43 

86.45 ± 3.68 

II 

Gaztanaga et al. (2009) 

0.44 

82.6 ±7.8 

II 

Blake et al. (2012) 

0.48 

97.0 ± 62.0 

I 

Stem et al. (2010) 

0.57 

92.4 ± 4.5 

II 

Samushia et al. (2013) 

0.593 

104.0 ± 13.0 

I 

Moresco et al. (2012) 

0.6 

87.9 ±6.1 

II 

Blake et al. (2012) 

0.68 

92.0 ± 8.0 

I 

Moresco et al. (2012) 

0.73 

97.3 ± 7.0 

II 

Blake et al. (2012) 

0.781 

105.0 ± 12.0 

I 

Moresco et al. (2012) 

0.875 

125.0 ± 17.0 

I 

Moresco et al. (2012) 

0.88 

90.0 ± 40.0 

I 

Stem et al. (2010) 

0.9 

117.0 ±23.0 

I 

Simon et al. (2005) 

1.037 

154.0 ±20.0 

I 

Moresco et al. (2012) 

1.3 

168.0 ± 17.0 

I 

Simon et al. (2005) 

1.363 

160.0 ±33.6 

I 

Moresco (2015) 

1.43 

177.0 ± 18.0 

I 

Simon et al. (2005) 

1.53 

140.0 ± 14.0 

I 

Simon et al. (2005) 

1.75 

202.0 ± 40.0 

I 

Simon et al. (2005) 

1.965 

186.5 ±50.4 

I 

Moresco (2015) 

2.3 

224.0 ± 8.0 

II 

Busca et al. (2013) 

2.34 

222.0 ± 7.0 

II 

Delubac etal. (2015) 

2.36 

226.0 ± 8.0 

II 

Font-Ribera et al. (2014) 


Note. — Here the unit of H{z) is km s~^ Mpc“^. “I” quoted 
in this table means that the H(z) value is deduced from the dif¬ 
ferential age method, whereas “H” corresponds to that obtained 
from the radial BAO method. 

used to constrain DE EOS at high redshifts. Simon et al. 
(2005) derived constraints on DE potential using H(z) results 
and supernova data. Samushia & Ratra (2006) began applying 
these measurements to constraining cosmological parameters 
in various DE models. In the meanwhile, DE evolution came 
into its own as an active research field in the last twenty years 
(Sahni & Starobinsky 2000; Carroll 2001; Peebles & Ratra 
2003; Copeland et al. 2006; Clifton et al. 2012). To sum up, 
the OHD are proved to be very promising towards understand¬ 
ing the nature of DE. 

In the next section, the OHD in Table 1 are used to 
test ACDM models under two-point Omh^(z2',Z\) diagnos¬ 
tic. Then, according to this diagnostic, we derive model- 
dependent values of Ho assuming that ACDM is valid. 

3. THE APPLICATIONS OE OMH^{Zr,Zi ) DIAGNOSTIC TO OHD 
3.1. The test of ACDM 

In our analysis, the validity of Om(z) diagnostic can be 
tested using H{z) results from cosmological independent mea¬ 
surements. Unlike the cases where only some special H{z) 
data from the radial BAO method are employed to test DE 
models (Sahni et al. 2014; Hu et al. 2014), thirty-six H(z) data 
points (as shown in Table 1) are considered in this work. Eor 
the sake of mitigating the contaminations from systematics 
of different H(z) observational methods performed by several 
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Fig. 1.— The full OHD set, binned data and standard flat ACDM model. 

The full and binned OHD are represented by gray and color-coded points, 
respectively. The yello\v cui've shows theoretical H{z) evolution based upon 
the adopted ACDM model which is described in Section 2. 

research groups, we bin all the OHD into four data points in 
four redshift ranges: 0.0 - 0.5, 0.5 - 1.0, 1.0 - 2.0, and > 2.0. 
We take an inverse variance weighted average of all the se¬ 
lected data in each redshift range in the following manner. 
Assuming that //,(z) represents the /th observational Hubble 
parameter data point with crHi(z) denoting its reported obser¬ 
vational uncertainty, in light of conventional data reduction 
techniques by Bevington & Robinson (2003, Chap. 4), it is 


straightforward to obtain 



Hiz)^ 

i:,(/7/(z)/cr^,,(,)) 

(4) 


4(.) - 

1 

(5) 



where Hiz) stands for the weighted mean Hubble parameter 
in the corresponding redshift range, and cr^(^) is its uncer¬ 
tainty. As a consequence, four weighted means of H{z) and 
corresponding uncertainties are computed to be 80.28 ±1.51 
kms-‘ Mpc-', 94.59 ± 2.76 kms-‘ Mpc^', 159.99 ± 7.89 
km s“' Mpc^', 223.81 ± 4.4 km s“' Mpc^*, at redshift zi - 
0.27, Z 2 - 0.73, Z 3 = 1.48, Z 4 - 2.33 respectively. These 
results are also shown in Figure 1. 

If Om{z2',Z\) is always a constant at any redshifts, then it 
demonstrates that the DE is of the cosmological constant na¬ 
ture. In order to compare directly with the results from CMB, 
Sahni et al. (2014) introduced a more convenient expression 
of the two-point diagnostic, i.e.. 


Omh^(z2',Zi) 


h^(z2) - h^(zi) 

(1 +Z2)^ -(1 +Zl)^’ 


(6) 


where h(z) = H(z)/IQQ km s“^ Mpc^'. The four H(z) points 
calculated based upon the aforementioned binning method 
therefore yields six model-independent measurements of the 
Omh^{z2',Zi) diagnostic, which are shown in Figure 2. 

In ACDM, we have Omh^ = The value of is 

constrained tightly by the Planck observations to be centered 
around 0.14 for the base ACDM model fit(Planck Collabora¬ 
tion et al. 2014): the Planck temperature power spectrum data 
alone gives 0.1423 + 0.0029, the Planck temperature data with 


Fig. 2.— The six Omh^ values computed from the four binned H(z) data. 
Different colors refer to different pairs of redshift bins. The uncertainty is 
estimated conforming to errors of each H{z) pair. The black line and shaded 
region con'espond to the best-fit value of and its uncertainties retrieved 
from the Planck CMB data. 

Planck lensing gives 0.1414 ± 0.0029, and the Planck temper¬ 
ature data with WMAP polarization at low multipoles gives 
0.1426 + 0.0025, all at 1 cr confidence level (CL). In this paper, 
we choose 0.1426 + 0.0025 as the Planck value. As shown in 
Figure 2, the testing pairs Omh^(zi', zj) and Omh^{z 2 ', Z 3 ) are 
consistent with ACDM at Icr CL, whereas for Omh^(z 2 ', Z 4 ) 
and Omh^izy, Z 4 ), ACDM can only be accommodated at 2cr 
CL. Particularly, for Omh^{zi,Z 2 ) and Omh^{zi',Z 4 ), ACDM 
is not compatible even at 2cr CL. Although deviations from 
ACDM exist for some pairs, cautiously, we cannot rule out 
the validity of ACDM. 


3.2. Measurements of Hq 

Recently, Hq determination remains in disagreement about 
the range of possible values. Some measurements of Hq are 
based on the calibrations from local distance indicators (Riess 
et al. 2011, 2012; Freedman et al. 2012; Chavez et al. 2012; 
Tammann & Reindl 2013; Fiorentino et al. 2013; Efstathiou 
2014). Others come from the global parameter fit to CMB 
anisotropy observations (Planck Collaboration et al. 2014) 
and extrapolation using existing Hiz) data (Busti et al. 2014). 
Some of the these results are also summarized in Figure 3. 

Unlike the previous subsection, here we derive Hq con¬ 
straints utilizing the two-point Omh^iz2',Z\) diagnostic, under 
the assumption that the values of Omh^iz 2 \Z\) are the same 
for every Hiz) pair, namely, ACDM being considered. The 
expression of the Omh^iz 2 ', Zi) diagnostic is thus modified into 


Omh^iz; 0) 


(7) 


A noteworthy characteristic of these Omh^ diagnostic is that 
it shall be a constant in ACDM. In order to determine this con¬ 
stant, similar to what we have done in the previous subsection, 
we take an inverse variance weighted average of the six Omh^ 
values from the four binned Hiz) points to get the final best-fit 
value and its uncertainty, i.e., Omh^ - 0.125 ± 0.004. Com¬ 
bining with Equation 7, it is then straightforward to derive the 
value of Hq (method (I)). In the meanwhile, we perform the 
same process except that the Omh^ value originates from the 
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Fig. 3.— DiflFerent determinations of Hq from different groups calculated 
using local calibrations or global parameter fits. The first ten results, from 
left to right, represent the studies of Riess et al. (2011), Riess et al. (2012), 
Freedman et al. (2012), Chavez et al. (2012), Tammann & Reindl (2013), 
Fiorentino et al. (2013), Efstathiou (2014) (with one distance anchor), Efs- 
tathioLi (2014) (with three distance anchors), Busti et al. (2014), and Planck 
Collaboration et al. (2014), respectively. The last two values are the results 
from method (I) and (II) presented in this work. 
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Fig. 4.— Hq values from the binned and observational H(z) data. Magenta 
points represent those derived from the inverse variance weighted method; 
meanwhile, black stars mark the results from the Planck measurement. The 
values of Hq, from left to right, are obtained from the thirty-six H{z) data 
points in ascending order of redshift. The insert in this figure shows derived 
Hq values from binned H{z) data with the same color coding. 


Planck measurement (method (II)). The results given by these 
two approaches are shown in Figure 4 for comparison. Gener¬ 
ally speaking, the values of Hq from the Planck measurements 
are lower than those from our weighted method. 

Utilizing methods (I) and (II) described above, we are able 
to determine all corresponding Hq values from each H(z) 
point (see Figure 4). However, in method (II), there are four 
observational H(z) data points which cannot be used to attain 
Ho value, due to the minus values of (^f )^ - Omh^{z', 0)((1 -H 
z)^ - 1). The associated redshifts of them are z - 1.53, 
z = 1.965, z = 2.34, and z = 2.36. We stmggle with the 
same problem when evaluating Hq from the fourth binned 
H{z) point. One cause may be ascribed to systematics induced 


et al. (2015) 

from observations and reduction procedures of both the H(z) 
data and the Planck measurement; there are probably other 
reasons latent in deep explanations, such as the deviation of 
ACDM. 

Then, we use the thirty-six Hq measurements (from method 
(I)) and thirty-two Hq measurements (from method (II)) to 
get the final weighted Hq values. Figure 3 summarizes the Hq 
results reported by different research groups. We obtain Hq - 
71.23 + 1.54 km s“^ Mpc^' in method (I) and Hq — 69.37 + 
1.59 km s“^ Mpc^' in method (II), both at Icr CL (see the last 
two results in Figure 3). Note that here the Hq values from 
the two methods are derived assuming an underlying ACDM 
model. 


4. UNCERTAINTY ON DE RECONSTRUCTION FROM ERROR OF 

OHD 

One of the key issues in cosmology is that whether the DE 
evolving its EOS or not, which reveals the underlying prop¬ 
erties of DE. The test of the deviation of the observational 
data from ACDM , as well as w(z) evolution with respect to 
redshift is difficult. One way to overcome this problem is to 
directly reconstruct the effective DE EOS via observations, in 
which errors need to be considered attentively. The section 
above indicated a way to test ACDM with OHD. In this sec¬ 
tion, we try to explore how the error in OHD propagate into 
w(z) at certain redshift during the reconstruction of w(z). 

A flat cosmological model with a DE term Qa = 1 - 
is implemented to analyze the error propagation. We first dis¬ 
card the errors in parameters but H{z) for simplicity. In this 
case, w(z) and its uncertainty are reformulated as 


“<’)=-A- 

l0g(l+8)(l ^ 

3 


( 8 ) 


^ ^ ^ 2H{z)6Hiz) 

6 w(z) = --r- 

3H2(^-OJl+z)3)ln(l+z) 


(9) 


Equations 8 and 9 indicate the relation between 6 H{z) and 
6 w(z). Here we fixed the relative errors of w(z) to be 20% 
and 50% arbitrarily to derive the required dH(z) to compare 
with the observed H(z) errors. The results are shown in fig¬ 
ure 5, where Hq - 67.1 km s“' Mpc^' and Q,„ = 0.3175 
are taken from the latest Planck constrains (Planck Collab¬ 
oration et al. 2014). It is worth noticing that high accuracy 
OHD is especially necessary in w(z) reconstruction at low 
redshift. As an example, the error of H{z) at z = 0.0708 
which is the nearest available OHD with the best fit value 
69.0 km s ' Mpc ', should be less than ±1.0 km s * Mpc ' 
and ±2.5 km s“' Mpc“* if we expect a reconstruction of w(z) 
with accuracy within 20% and 50% respectively. It should 
be noted, though, that there are six OHD points which cannot 
be applied in w(z) reconstruction, due to the minus values of 

- Qm(l + z)^ (that means > 1). The respective red- 

shifts of them are z = 0.88, z = 1.53, z = 1.965, z = 2.3, 
z = 2.34, and z = 2.36. 

However, the absence of errors in other parameters is of 
course an idealization. We then consider the more complex 
situation that the uncertainty of w(z) not only come from the 
error of H(z), but also have a bearing on the errors of Hq and 
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Fig. 5.— The ratio of required H{z) error to the best fit value of OHD when 
the relative error of vy(z) is fixed at 20% or 50% with ideal assumption that 
the uncertainty of reconstructing w{z) only comes from the prior on the error 
of H(z). The gray error bar represents the ratio of observational error of H{z) 
data to its best fit value. Red error bar means the ratio of predicted error 
of H{z) data to its best fit value if 6w{z)lw{z) is fixed at 20%. Cyan color 
represents the situation of 6w{z)lw{z) = 50%. 

Q,„. What deserves special mention is that we do not con¬ 
sider the correlation between parameters owing to the lack of 
necessary information. In this case, the uncertainty of w(z) is 
given by 


6 w(z) = 


2 H(z) 


3//2(^-QJl+z)3)in(l+z) 


-2H\z)H^ 


-3 


3(^-n„(l+z)3)ln(l+z) 
_ (1+z)^ _ 

3(^-n„(l+z)3)ln(l+z) 


6^mz) 


6^Ho 


( 10 ) 


1 


3(l-n„,)ln(l+z) 




Here we choose Hq - 67.4 +1.4 km s“* Mpc^' and Qm = 
0.314 + 0.02 as the prior which are constrained tightly by the 
Planck observation (Planck Collaboration et al. 2014). Ob¬ 
viously, increased accuracy OHD is required if more param¬ 
eters are considered to contribute to the uncertainty of w(z). 
Reconsider the case of H(z) at z=0.0708, when the relative 
uncertainty of w(z) is constrained at 50%, we could derive 
the required 6 H(z) - ±2.2 km s“' Mpc^', however, required 
6 H(z) cannot be derived when 6 w(z)/w(z) is fixed at 20% be¬ 
cause the relative uncertainty of w{z) which comes from the 
errors of Hq and is already over 20%. Some data points 
in figure 6 do not have a required error bar (at z - 0.0708, 
0.09, 0.17, and 0.28 for dw(z)/w(z) = 20%; z = 0.17 for 
6 w{z)lw{z) - 50%) comparing with that of in figure 5 due to 
the same reason. 


5. CONCLUSIONS AND DISCUSSIONS 

In this paper, motivated by the investigation on the nature 
of DE, we have tested the validity of ACDM by the two-point 


Fig. 6.— Same as Figure 4, except that the uncertainty on the reconstructed 
w(z) come from the priors on the errors of H{z), Hq, and Magenta er¬ 
ror bar means the ratio of predicted error of H{z) to its best fit value if the 
uncertainty of w(z) is fixed at 20%. Blue color represents the situation of 
50%. 

Omh^{z 2 ',Zi) diagnostic using thirty-six OHD which are ob¬ 
tained from differential age and BAO methods. In order to 
obtain more reliable result, we binned all the OHD into four 
data points according to the redshift distribution to mitigate 
the observational contaminations. Moreover, in order to de¬ 
termine the best fit value and the uncertainty of binned Hiz), 
we took an inverse variance weighted average method. Our 
result furthers a discussion about the essence of DE, namely, 
we cannot rule out the validity of ACDM relying on nothing 
more than the current OHD. 

In a natural way, on the premise of ACDM, we further 
used the two-point Omh^izi',Zi) diagnostic to derive the Hub¬ 
ble parameter at z = 0. Eor more comprehensive analysis, 
we adopted two methods to determine the value of Omh^, 
one derived from binned H(z) points (method (I)), the other 
was based on the result obtained from Planck measurement 
(method (II)). Throughout above two methods, corresponding 
Hq value has been reconstructed from each observational H(z) 
point logically. However, there are four OHD at z = 1.53, 
z = 1.965, z = 2.34, z = 2.36 cannot be used to attain the 
Hq value by method (II). One reason may be attribute to sys- 
tematics come from observations and reduction procedures; 
there can be others deep reasons, for example, the deviation 
of ACDM. Thus, attention should be paid to some thorough 
effort on improving spectra quality to yield more H(z) deter¬ 
minations with small error in future work. 

Utilizing the Hq determinations with error corresponding 
to each H{z) data derived by method (I) and (II), we reached 
the final weighted Hq - 11.23 ± 1.54 km s“' Mpc^' (method 
(I)) and Ho = 69.37 + 1.59 km s^' Mpc^^ (method (II)). Ad¬ 
mittedly, our results about Hq is tentative as well as model- 
dependent. The conclusion is obtained under the assumption 
of ACDM model. Therefore, care must be taken not extrap¬ 
olating our conclusions beyond this assumption, especially 
when used in constraining cosmological parameters. 

Eurthermore, in the end of this paper we explored how the 
error of H{z) impact on the uncertainty of w(z). To illustrate 
the effect from H(z) data, it is worth making a comparison be¬ 
tween the two assumptions, namely, the uncertainty on w(z) 
simply comes from the error of H(z), or also involves the con¬ 
tributions from other parameters’s uncertainties. We argue 
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that the current precision on OHD is not sufficient small to 
ensure the reconstruction of w(z) in an acceptable error range, 
especially at the low redshift. In other words, the ability of 
available individual OHD on the DE evolutionary study is 
very limited. 

Finally, we note that it is necessary to place our main em¬ 
phasis of research on BAO observation. Using BAO peak po¬ 
sition as a standard ruler in the radial direction or combin¬ 
ing measurements of BAO peak and Alcock-Paczynski dis¬ 
tortion, we can limit the precision of H{z) data better than 
7% (Gaztanaga et al. 2009; Blake et al. 2012). As putting 
into operation of future space and ground-based telescopes 
(James Webb Space Telescope, Wide-Field Infrared Survey 
Telescope, planned adaptive optics systems with Keck, Farge 
Synoptic Survey Telescope, and Thirty Meter Telescope et 
al.), more high-redshift, high-accuracy H{z) determinations 
from BAO observations will undoubtedly perform a very use¬ 
ful role in the future study of the DE. 
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